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ABSTRACT 

We use new Suzaku observations of PKS 0745-191 to measure the thermodynamic properties 
of its ICM out to and beyond I200 (reaching 1.25r2oo) with better accuracy than previously 
achieved, owing to a more accurate and better understood background model. We investigate 
and resolve the tensions between the previous Suzaku and ROSAT results for PKS 0745-191, 
which are found to be principally caused by incorrect background modelling in the previous 
Suzaku analysis. We investigate in depth the systematic errors affecting this observation, and 
present temperature, density, entropy and gas mass fraction profiles reaching out to and be- 
yond the virial radius. We find that the entropy profile flattens in the outskirts as originally 
observed in the previous Suzaku analysis, but that the flattening starts at larger radius. The 
flattening of the entropy profile and our mass analysis suggests that outside ~17' (~1.9 Mpc) 
the ICM is out of hydrostatic equilibrium or the presence of significant non-thermal pressure 
support. 
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1 INTRODUCTION 

Accurate measurements of the intracluster medium (ICM) of 
galaxy clusters to the virial radius (taken to be r2ocEl)! impor- 
tant for many reasons. The virial radius represents the boundary 
within which the cluster ICM is expected to be in hydrostatic equi- 
librium, and outside of which matter is still accreting onto the clus- 
ter as it continues to form. Studying cluster outskirts allows the 
cluster formation process to be better unde rstood, thus constrain- 
ing simulations of cluster formation (such as lRoncarelli et al.l2006l 
lBumsetal.l|201 (j) while allowing us to constrain models for the 
baryon fractions of clusters (for example lYoung et al.lioT l]) which 
are dependent on non-gravitational processes such as feedback and 
cooling. 

We can also investigate where the assumption of hydrostatic 
equilibrium breaks down. The hydrostatic equilibrium assumption 
is important because it is used when calculating the masses of 
galaxy clusters, as we see later in section|5] Accurate cluster masses 
are required for using galaxy clusters as probes of cosmological 
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r200 is the radius witliin which the mean density of the cluster is 200 
times the critical density required for a fiat universe, pc - The mass enclosed 
within r200 is M200=4/37rpcr|Q(,. 



parameters using the cluster mass function JVikhlinin et"al]l2009t) . 
and understanding the gas mass fraction to the virial radius im- 
proves the_reliability_with which it can be used as a cosmological 
probe jAllen et al.l2008[) . 

In lGeorge etZI ( l2009h (hereafter G09) the ICM of the galaxy 
cluster PKS 0745-191 was investigated to beyond the virial radius 
they calculated for the first time. This exploited Suzaku' s low in- 
strumental background, resulting from its low orbit which places it 
inside the Earth's magnetopause, allowing the low surface bright- 
ness outskirts to be probed. This development was made possible 
by the fact that PKS 0745-191 is very X-ray bright (it is the bright- 
est cluster beyond z = 0. 1). 

One of the most interesting findings of G09 was that the 
entropy profile of the ICM flattened in the outskirts outside 
10'=l.lMpc, deviating significantly from the theoretical prediction 
of a r^'^ power law incr ease assuming gravitational structure for- 
mation jVoit et alJ|200^ . However the location of PKS 0745-191 
just 3 degrees above the galactic midplane means that it has a com- 
plicated and spatially variable background and absorbing column 
density which could not be properly understood due to the lack of 
suitable background pointings when G09 was written. In order to 
fully understand the properties of PKS 0745-191 in the outskirts, an 
accurate background model is needed whose errors are thoroughly 
understood. 
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Further clusters have since b een studied to the v irial ra- 
dius with Suzakij [Abe ll 2204 l lReiprich et al] |2009|) ; Abe ll 
1975 teautz et alJ l2009h: Abel l 1413 dHoshino et alj l2010h : 
Abell 1689 iKawahara da et alj I2OI0I) ; the Perseus cluster 
dSimionesciTetal. (2011a)); Ab ell 2142 jAkamatsu et al l I2OIII) : 
Abell 2029 l l Walker et alJT2012h and Hydra A dSato et al.ll2012h ]. 
and a statistical sample with general trends has started to emerge. 
One common feature that is observed for most of these clusters 
is the flattening of the entropy profile from around 0.5r2oo out- 
w ards, causing the entropy to be below the powerlaw prediction 
of IVoit et alj l l2005l) in the outskirts. One explanation of this is 
that the cluster ICM is out of hydrostatic equilibrium in the out- 
skirts, and that non-thermal pressure support from turbulence and 
bu lk motions is impo rtant, as is found in the numerical simulations 
of iBums et al However the results for the Perseus clus- 

ter (which because of its closeness and thus large angular extent 
has been studied with better resolution than the other clusters) sug- 
gest that the entropy profile flattening is a consequence of the gas 
density in the outskirts being overestimated due to gas clumping, 
which also causes the measurements of the gas mass fraction and 
gas pressure in the outskirts to be higher than expected. 

Numerical simulations of galaxy cluster outskirts indicate that 
the contributions of non-thermal pressure support increase with ra- 
dius in relaxed clusters, leading to an underestimate of the hydro- 
static mass if only the thermal pressure (which is measured directly 
in X-ray observations) is used in the hydrostatic equilibrium equa- 
tion. Random gas motions become increasing important in the out- 
skirts of galaxy cluster (from rsoo outwards) due to merging ac- 
tivity, the continual accretion of gas onto the cluster along large 
scale structure filaments, and due to su personic motions of galaxies 
through the ICM. The simulations of iLau et al. I ( l2009h suggested 
that the non-theiTnal pressure contribution in relaxed clusters in- 
creases with radius to around 20-30 percent at r2oo, resulting in the 
hydrostatic mass derived from the thermal pressure underestimat- 
ing the actual mass by '--^lO percent at rsoo, rising to ~ 20 percent 
at r2oo. Similar results for the mass und erestimate in the c luster 
outskirts were found in the simulations of iBums et al. I (|20T1 . who 
found that at the thermal energy density decreases with radius while 
the energy density in turbulent kinetic energy increases with radius, 
with both contributions being approximately equal at ~40 percent 
each at r2oo, with the remaining ^--^20 percent being provided by 
radial kinetic energy due to infalling gas. 

Simulations also indicate that gas clumping plays an impor- 
tant role in effecting t he measured gas d ensity in the outskirts 
jMathiesen erZlll999l iNagai & Laulboilh . This occurs because 
the ICM emits X-rays due to thermal bremsstrahlung, with the flux 
being proportional n|f, so if the ICM is clumped then this leads 
to more X-ray emission than if it is uniform with the same aver- 
age density because (n|f) > (n£f)^. Therefore, if the ICM is as- 
sumed to be unform when it is in reality clumped, the gas density 
will be overestimated, leading to an overestimate of the gas mass 
fraction and the gas pressure, and an underestimate of the gas en- 
tropy (S=kT/ne''''). In lNagai & Lad 1 I2OIII) gas clumping is found 
to be most important outside r2oo, though the scatter between dif- 
ferent simulated clusters is large. The exact nature of the clumping 
(its size scale and distribution, its azimuthal variation and relation 
to large scale structure filaments) remains uncertain both in simu- 
lations (where it depends sensitively on the precise ICM physics, 
such as the rate of cooling and star formation) and in observation. 

Tensions between the Siizaku results of G09 and mea- 
surements wit h other X-ray obse rvatories have recently been 
highlighted in lEckert et al. (201 lall (hereafter Ella), which used 



ROSAT PSPC observations, and found density and surface bright- 
ness profiles which are inconsistent with those presented in G09 
in the outskirts. In addition the temperature profile found in G09 
is much lower than that o btained in the 700-1500 kpc reg ion us- 
ing X MM-Newton CSnowden et al.ll2008h and Swift (More tti et al.1 
|20I1|) . and the mass determination in G09 is lower than and incon- 
sistent with previous measur ements (see table |4]for a comparison). 
In addition. lSato et al] ^20121) has shown that while the entropy pro- 
files (scaled by the average ICM temperature) for the clusters inves- 
tigated to r2oo with Suzaku so far are consistent with one other with 
low scatter, the entropy profile found in G09 for PKS 0745-191 is 
lower in the outskirts and inconsistent with the other profiles. 

Ella attributed these tensions to incorrect background mod- 
elling. As no nearby background pointings using Suzaku were 
available when G09 was written, the background was modelled us- 
ing data from the Lockman Hole, which was found from Rosat All 
Sky Survey (RASS) data to have a similar count rate in the 1.0-2.0 
keV band as the regions surrounding PKS 0745-191. The Lockman 
Hole data were also chosen because they were taken at a similar 
time as the original PKS 0745-191 data, thus minimising the ef- 
fects of time variability of the contamination on the optical block- 
ing filter (OBF) in front of the detectors (which was not completely 
understood when G09 was written). However the Lockman Hole 
is located at a much higher galactic latitude than PKS 0745-191, 
meaning that the galactic column density is much lower, and that 
the source of the 1-2 keV emission observed in the Lockman Hole 
data would not be expected to penetrate the absorbing column at the 
latitude of PKS 0745-191. Another contribution to the background 
at low galactic latitude must exist to compensate for the decrease 
in flux resulting from greater absorption, and its spectrum may be 
significantly different. 

To further study the properties of the cluster outskirts and to 
better constrain our model of the background flux, we have ob- 
tained new Suzaku observations near PKS 0745-191. In this paper 
we analyse these new data along with the original data used in G09 
and other archival observations to investigate the source of incon- 
sistency seen in previous work. 



2 OBSERVATIONS AND DATA REDUCTION 

To supplement the original 5 pointings in the cross formation used 
in G09, which were taken in 2007 May, five more pointings were 
taken in 2010 May and 2010 November, which extend the cross 
to the North, West and South and fill in the comers between these 
three arms, as shown in Fig.[T] This allows a local background de- 
termination to be made and allows spatial variations in the back- 
ground parameters to be better understood, which is important for 
quantifying the uncertainty in the background and its affect on the 
cluster properties derived from spectral fits. 

The data were reduced using the same method as that de- 
scribed in lWalker et al. 1( |2012|) . The ten Suzaku XIS (X -ray Imaging 
Spectrometer) pointings of PKS 0745-191 were reduced using the 
ftool AEPIPELINE which perfoiTns the standard filtering described 
in section 6.2 of the Suzaku ABC guidfQ using the latest calibra- 
tions in HEASoft version 6.11.1. The screening criterion COR>6 
(which removes time intervals where the cut off rigidity was below 
6 GeV) was included which reduces the good time interval (GTI) 
total exposure time to the values shown in Table [T] We generated 
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Figure 1. Le/fiMosaic of Suzaku pointings. Pointings used in G09 are in the central cross, and tlie new pointings extending tlie cross are sliown in red. Annuli 
are between 0.0'-2.5', 2.5'-6.0', 6.0'-8.8', 8.8'-11.7', 11.7'-14.5', 14.5'-17.3', 17.3'-20.2', 20.2'-23', 23'-28', 28'-32' and 32'-41'. Right: Point sources 
removed from tlie analysis using a threshold flux of 10^^'^ erg cm^^ s^^ in the 2-10 keV band. In both panels the vignetting has not been corrected for, and 
the images are for illustrative purposes only. 
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Figure 2. Checking for SWCX contamination during the periods of observat ion using the 0+ '^/0 +^ ratio obtained with the ACE spacecraft. Contamination is 
important when O'^^/O'^^ is greater then ~ 0.2. Left: Original data used in jOeorge et al.ll2009l) Centre: Pointings taken in May 2010. Right: Pointings taken 
in November 2010. 



exposure maps using XISEXPMAPGEN to identify the regions near 
the edge of the detectors which are removed due to their low effec- 
tive area. The calibration regions at the comer of the XIS detectors 
were removed. 

Non X-ray background (NXB) spectra were created using 
XISNXBGEN and were subtracted from the spectra when perform- 
ing fits. Redistribution matrix files (RMFs) were made using XIS- 
RMFGEN and Ancilliary Response Files (ARFs) were made using 
XISSIMARFGEN. We use two ARFs: one assuming a uniform circu- 
lar source of radius 20' (larger than the detector field of view, FOV) 
which is used for the uniform background model, and one made us- 
ing the j3 model best fit to the surface brightness profile f ound with 
our Suzaku data using the method o f iBautz et alj ( |2009|) . which is 
used for the cluster emission model. 

We remove point sources identified by eye in the Suzaku im- 
ages as shown in Fig. [T] meaning that we resolve the cosmic X- 



ray background (CXB) as point sources to the same threshold flux 
(1 X 10~^^ ergs cm~^ s~^) across all of the poi ntin gs (as is done 
in 'B autz et all l2009l Isimionescu et al. (201 lail and IWalker et all 
2012). Point sources below this flux are unresolved and are mod- 
elled as a powerlaw component of the background. We derive the 
expected unresolved CXB flux and its variance for the annular ex- 
traction regions later in section |33] 

All of the light curves were checked to ensure that no flaring 
had occurr ed during the observati on in the energy band studied. 
Following iHumphrev et al.l l l201lb we obtained ACE (Advanced 
Composition Explorer) data using the SWICS (Solar Wind Ion 
Composition Spectrometer) instrument for the O'^'^/O'^^ ratio at 
the times of the observations to gauge the effect of solar wind 
charge exch ange (SWCX) on our observations, following the ob- 
servation of ISnowden et al.l ( |2004|) that SWCX is negligible when 
Q+7/Q+6 jg^^ jjjgjj ^ Q 3 shown in Fig. |2] we find that 
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O^'^/O^® never increases significantly above ~ 0.2 during the ob- 
servations, and conclude that SWCX contamination is negligible. 
As a further check we filter out the time intervals where O^^/O'*'® 
is greater than 0.2 and find that the spectral fits to the data are not 
affected. 

All fits were performed in XSPEC 12.7.1b using the extended 
C-statistic. We only use data from the front illuminated XIS de- 
tectors (XISO and XIS3) due to the higher NXB level of the back 
illuminated XISl detector. We fit the XISO and XIS 3 data for both 
of the editing modes (3x3 and 5x5) simultaneously (so each fit is 
done to 4 spectra simultaneously). 



3 BACKGROUND SUBTRACTION AND MODELLING 

An accurate understanding of the outskirts of galaxy clusters re- 
quires the level, spectral shape and spatial variation of the X-ray 
background to be understood. The extragalactic component of the 
X-ray background (the cosmic X-ray background, CXB) which 
originates from unresolved point sources, and the soft foreground 
emission from the galaxy both need to be modelled accurately and 
the magnitude of their expected spatial variations over the cluster 
properly quantified. 

We treat each of these background components separately 
later, but first describe the problems with the background modelling 
used in G09 and describe the complicated absorbing hydrogen col- 
umn density in the area around PKS 0745-191 which has not pre- 
viously been taken into account. 



3.1 Previous Background Subtraction and Modelling 

In G09 the Lockman Hole was modelled as consisting of a single 
unabsorbed thermal component (a MEKAL model) at kT=0. 1 keV 
representing soft galactic emission, added to an absorbed power- 
law of index 1.4 representing the cosmic X-ray background (CXB) 
from unresolved point sources. This was motivated by the similar 
count rates of the Lockman Hole and the background region around 
PKS 0745-191 in the 1-2 keV band as determined from ROSAT All 
Sky Survey (RASS) data. However as PKS 0745-191 lies at a much 
lower galactic latitude than the Lockman Hole [for PKS 0745-191 
(l,b)=(236.4,3.0) and for the Lockman Hole (l,b)=(149.7,53.2)], the 
absorbing column density is higher in the direction of PKS 0745- 
191 (with an average of 4.2 x 10^^ cm^, compared to 6.2x10^'' 
cm^ for the Lockman Hole). This means that the actual galactic 
emission must be higher than that in the direction of the Lock- 
man Hole (to compensate for the higher column density and yield 
the same count rate), which may be the result of a component not 
present in the Lockman Hole direction. 

In Fig. [3] we compare the actual background spectrum from 
the background regions in our new Suzaku observations (regions 
outside 23') with the Lockman Hole background level used in G09. 
The CXB level used in G09 (7.26 x 10"* photons cm"^ s"^ at 1 
keV for a circular area of radius 20 arcmins) corresponds to a 2-10 
keV band flux of 1.37x 10"^^ ergs cm"^ s"^, which as we find 
in section [33] is at the lower limit of what we would expect given 
a constant threshold flux of removed point sources of 10~^^ ergs 
cm~^ s~^, causing this background component to be systemati- 
cally underestimated. This is likely because G09 did not use a con- 
stant threshold flux for resolving point sources, and because only 
one background pointing was used whose area was insufficient to 
cover a representative sample of point sources in the CXB. As the 
fit was only performed in the 1-5 keV band, the fit below I keV is 
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Figure 3. Comparison of the Lockman Hole background model used in 
G09 with the background spectral data points from the Suzaku background 
regions (23' to 41') of our PKS 0745-191 pointings. The causes of the dis- 
crepancies between the model and the data are described in the main text. 



poor indicating that the background model for the soft emission is 
inco rrect. 

iMasui et alj ^2009*) examined a background pointing in the 
galactic midplane, (l,b)=(235. 00,0.00), 3 degrees away from PKS 
0745-191, and deduced a significantly different backg round model 
to G09 resulting from its lower galactic latitude. In iMasui et"ail 
(2009) three background models are presented using two different 
hypotheses for the background which attempt to explain a narrow 
bump in the spectrum which peaks at ~ 0.9 keV. The first type, 
(which is explored in their model A and model B) uses the same 
CXB and galactic local hot bubble (LHB) model as G09, but in- 
cludes higher temperature thermal components (for model A one 
VAPEC component at 0.766keV is used, while in model B a VAPEC 
component at 0.658 keV is added to an APEC component at 1.5 
keV). In the preferred model (their model C), the hot gas com- 
ponent is replaced with an absorbed bremsstrahlung component at 
kTc=0. 183 keV which they use to model the contribution of un- 
resolved faint dM stars to the spectrum (whose number density is 
higher close to the galactic plane), which is claimed to give a more 
physical model for the spectrum. 

From Fig.[3]we see no significant evidenc e for th e bump in the 
spectrum at 0.9keV observed in Masui et al] i2009|), and we find 
that the background can be best fit assuming a model consisting of 
two thermal components (APEC) each of solar metallicity and zero 
redshift at 0.1 keV (unabsorbed, modelling the LHB) and 0.6 keV 
(absorbed), and this fit is shown in Fig. lA2l This is most likely be- 
cause the dM star component which contributes this background 
com ponent decreases r apidly away from the galactic plane (see Fig 
4. of i Masui et id][2009l which shows the dM star contribution de- 
creasing by a factor of ~ 10 from its galactic midplane value at 
the latitude of PKS 0745-191, b=3.0), and we find no statistically 
significant evidence for its presence. 



3.2 Galactic hydrogen column density variations 

The column density of hydrogen in the direction of PKS 0745- 
191 is high due to its location near to (3 degrees above) the 
gala ctic plane, and from the Leiden/ Argentine/Bonn (LAB) sur- 
vey jKalberla etai]|2005l) and the Dickey and Lockman survey 
l lDickev & LockmanI [T990l) is known to be highly spatially vari- 
able, though due to the course sampling of these HI surveys the 
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Table 1. Observational parameters of the pointings 



Name 


Obs. ID 


Position 


Total exposure per detector (ks) 
(3x3 and 5x5 editing modes) 


RA 


Dec (J2000) 


10 


802062010 


Centre 


27.9 


116.8852 


-19.2901 


20 


802062020 


W 


26.0 


116.6543 


-19.2063 


30 


802062030 


N 


26.3 


116.9737 


-19.0727 


40 


802062040 


E 


30.6 


117.1155 


-19.3739 


50 


802062050 


S 


31.2 


116.7966 


-19.5079 


53 


805083010 


FarN 


31.8 


117.0925 


-18.8008 


54 


805084010 


NW 


19.3 


116.7033 


-18.9048 


55 


805085010 


FarW 


31.1 


116.3632 


-19.0871 


56 


805086010 


SW 


23.9 


116.4685 


-19.4632 


57 


805087010 


FarS 


32.4 


116.6522 


-19.7838 



full extent of the spatial variation is poorly understood. Infrared 
dust maps such as the 100 micron IRAS dust map shown in Fig.|4] 
which are known to correlate with column density variations, show 
a complicated environment not fully represented by the LAB and 
Dickey and Lockman survey values which are also shown in Fig. 
|4] The general trend around PKS 0745-191 shown in Fig.|4]is for 
the column density to increase from east to west, with a range from 
0.32 to 0.52 X 10^^ cm"'^ (using the LAB values), and the average 
value over the cluster is 0.42 x 10^^ cm~^. 

Previous X-ray observations of the bright centre of PKS 0745- 
191 have found column de nsities at the l ower end of the range in- 
dicated by the HI surveys. lAUerJ ( l2000h studied a ROSAT PSPC 
pointing of PKS 0745-191 and found a best fittin g column den- 
sity at the core of ~ 0.35 x 10^^ cm~^. IChen et alj lf2003:l studied 
XMM-Newton pointings and found a best fit central column density 
of 0.36 X 10^^ cm"^ with the EPIC-PN detector data and 0.4 x 
10^^ cm"^ with the EPIC-MOS data (this difference is attributed 
to a systematic calibration problem with these detector's response 
matrices). As we find later when we allow the column density to 
be a free parameter in our spectral fits, we find that the column 
density near the core is consistent with these results from ROSAT 
and XMM-Newton, but that away from the core the column density 
varies spatially (but is still consistent with the range from the HI 
radio surveys). This is shown in the right panel of Fig.|4] 

3.3 The CXB 

The CXB, which consists of unresolved point sources, has been 
studied extensively and progressively deeper observations have al- 
lowed the point sources fro m which it is compris ed to be resolved 
to progressively lower flux. iMoretti et ^ (l2003 h studied ROSAT, 
Chandra and XMM-Newton observations and obtained empirical 
best fits to the cumulative flux distribution (the number of sources 
have a flux greater than S, N(>S)) using a two powerlaw model. 
For studies of the outskirts of galaxy clusters, spatial variations 
in the CXB due to spatial variations in the number of unresolved 
point sources presents one factor limiting our ability to detect faint 
cluster emission in the outskirts. The deeper the CXB can be re- 
solved, the lower the number of unresolved point sources and thus 
the lower the spatial variations expected in the CXB due to spa- 
tial variations in the number of unresolved point sources, making it 
possible to measure fainter cluster emission with greater statistical 
significance. 

When modelling the expected CXB level and its expected 
level of fluctuation in the spectral extraction regions, we follow the 



approach o f lWalkeretal1 ( l2012h . For the total unresolved CXB flux 
in the 2-10 keV band, we use the value obtained in lMoretti et al.l 
( 2009) using Swift (2.18 ±0.13 x 10"" erg cm"^ s"^ deg'^), 
which is in complete agreement with the value determined in 
iMoretti et alj ( l20()3i) (2.24 ± 0.1 6 x 10~" erg cm~^ s ~^ deg"^). 
We avoid using ASCA results of iKushino etal] j2002l) due to the 
stray light problems from which ASCA suffered. 

To determine the amount of the CXB resolved, we integrate 
the cumulative flux distribution above the threshold flux to which 
our Suzaku observations allow point sources to be resolved (Sexci = 
10^^'' erg cm^^ s^^ in the 2-10 keV band) and subtract this from 
the total 2-10 keV band flux (in ergs cm~^ s~^ deg~^) as shown 
below, 

FcxB = 2.18 ± 0.13 X 10'" - j xSdS (1) 

which leads to an unresolved 2-10 keV band CXB flux of 
1.87±0.13 X I Q-" ergcm"^ s~ ^ deg-^ 

Following i Bautz et al.l j2009 l). we then calculate the expected 
spatial variations in the 2-10 keV band CXB flux due to the unre- 
solved point sources in each annular spectral extraction region (of 
solid angle Q.) using 

•jhxB = im) (^) X dS (2) 

from which we can see that the variance can be minimised by 
using larger extraction regions and a smaller threshold flux (Soxci)- 

The resulting one sigma uncertainties for each annulus in 
which we investigate cluster emission are tabulated in table [2] 
and these are the ranges by which we must vary the CXB back- 
ground model to understand the systematic error on the ICM 
measurements of our uncertainty on the CXB level. These pre- 
dicted variations are entirely consistent with thos e derived b y 
scaling the fluctuations observed by Ginga ( iHoshino et alj|201oh . 
and in the Suzaku backgro und observations of the Perseus cluster 
dSimionescu et al. (2011bll) . 

3.4 Modelling soft foreground. 

The soft foreground resulting from galactic emission is challenging 
to measure for PKS 0745-191 due to its location near the galactic 
plane and the highly spatially variable hydrogen column density in- 
dicated by the HI survey maps and the IRAS 100 micron dust maps. 
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Figure 4. Left: IRAS 100 micron map with LAB and Dickey and Lockman liydrogen column density values (units are 10^^ cm~^) overplotted (LAB results 
are first, then the Dickey and Lockman results). The black dashed circle is at 41' and represents the outer limit of all the pointings. The dashed blue circle of 
radius 23' is the same in both the left and right panels and represents the boundary beyond which we detect no statistically significant cluster emission. Right: 
Column density values (units are 10^^ cm^^) obtained by letting the column density be a free parameter in our Suzaku fits. The error on each value (combined 
systematic and statistical error) is ±0.05 X 10^^ cm~^. 



Table 2. Expected CXB one sigma vaiiations in each annulus in the 2-10 keV band (10 erg s ^ cm' -2 deg-2) 



0.0' -2.5' 


2.5'-6.0' 


6.0'-8.8 


' 8.8'-11.7' 


11.7'-14.5' 


14.5'-17.3' 


17.3'-20.2' 


20.2' -23' 


6.7 


3.1 


3.1 


2.5 


2.5 


2.4 


2.1 


2.3 



Therefore, in addition to spatial variation of the galactic emission, 
there is also spatial variation of the column density absorbing this 
emission. 

When fitting for the soft foreground using the Suzaku observa- 
tions outside 23', we find it is best modelled using two APEC com- 
ponents each with solar metallicity and zero redshift (one at 0.6 keV 
which is absorbed by the galactic column, and one at 0.1 keV which 
is unabsorbed and corresponds to emission from the Local Hot 
Bubble). The fit is not improved by including another low temper- 
ature absorbed component correspo nding to the hot halo surrou nd- 
ing the galaxy which is often used ISimionescu et al. (2011bll) . In 
general it is not possible to fit for the APEC normalisations of the 
galactic background components and the column density simulta- 
neously owing to the strong degeneracy between them, in that the 
same background emission level is achieved by reducing both the 
column density (i.e. less absorption) and the galactic APEC normal- 
isations (i.e. less galactic emission). 

To begin with we fit the 32'-4l' regions of the Suzaku point- 
ings with the background model described earlier, initially fixing 
the column density to the mean value over the cluster of 0.42 x 
10^^ cm~^ and allowing the APEC and CXB normalisations to be 
free parameters. We then compare the 23 '-32' regions in the new 
pointings with them to see whether we can see evidence for clus- 
ter emission. We find fluctuations in the soft band flux, but these 
can be explained by variations in the column density which are in 



the range expected from the LAB survey values. So there is no sta- 
tistically significant evidence for excess emission above the back- 
ground in the region outside 23'. 

Thus we find that although the soft background varies between 
23' and 41', it is consistent with the simple possibility of uniform 
galactic emission absorbed by a varying column density within the 
range of the LAB survey values. However the possibility remains 
that the intrinsic galactic emission and the column density are both 
varying, so for example it is not possible to say with certainty that 
a background region with higher than average soft band flux is a 
region with lower than average galactic column; it could instead 
be a region with higher than average galactic emission with the 
same absorbing column as other regions. It is not possible to break 
this degeneracy due to the courseness and uncertainties of the LAB 
survey galactic hydrogen column values. 

Due to this, we adopt a conservative approach when mod- 
elling the soft foreground emission which fully encompasses 
our uncertainty in the spatial variation of the galactic emis- 
sion and the spatial variation of the column density absorbing 
it. The Suzaku background regions are divided into annuli be- 
tween 23'-28', 28'-32' and 32'-4l' (yielding 13 regions in total 
as shown in Fig. [TJ so that the extraction regions are approx- 
imately the same size as the cluster emission annuli, allowing 
us to understand the expected variation of the soft foreground 
on this distance scale. We fit each Suzaku background region 
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with the phabs(apec(0.6keV)+powerlaw)+apec(0. IkeV) model us- 
ing the minimum absorbing column indicated by the LAB survey 
(0.32 X 10^^ cm~^) and then for the highest absorbing column 
(0.52 X 10^^ cm~'^). The galactic APEC and powerlaw noms are 
free parameters which we fit for. From this we obtain the maxi- 
mum and minimum galactic APEC norms for each component (the 
highest coming for the higher absorbing column), and this gives 
the full range within which we can say with certainty the intrin- 
sic galactic APEC norms must lie, and this range folds together our 
uncertainty in the spatial variation of the absorbing column and un- 
certainty in the intrinsic galactic emission. For all of the fits the 
CXB level is consistent within the expected va riance w i th that cal- 
culated in section [331 using the results of Mor etti et al.l l l2003h and 
the point source threshold flux of 10^^"^ erg cm~^ s"'^ in the 2-10 
keV band. 

As a further check we perform spectral fits to the 25 '-50' an- 
nulus of the ROSAT PSPC pointing for PKS 0745-191 (obs ID 
rp800623n00, which is shown in Fig. IA3t . dividing it into 4 sec- 
tors to the north, south, east and west, allowing us to measure the 
soft background in the east which is not covered by our Suzaku 
background pointings. The ROSAT data were reduced and filtered 
using the standard procedures described in the document 'ROSAT 
data analysis using xselect and ftools' 0. Particle background spec- 
tra were modelled using the ftool PCPARPHA and subtracted, and 
ROSAT ARFs were created using the ftool PCARF. The galactic 
APEC components are found to be in agreement with the values 
from the Suzaku background pointings, and the background level 
of the east direction is consistent with that of the other directions. 

We take the best fitting galactic APEC norms to be the values 
obtained when fitting all of the 13 background regions simultane- 
ously with a column density of 0.42x 10^^ cm^^ (the mean column 
density from the LAB survey) and these are norm(0.1 keV)=0.015 
cm^"" and norm(0.6 keV)=0.0013 cm^'', where the normalisations 
have units of 10""(47r)"^D;^^(l + z)'^ J nenudV , where Da 
is the angular size distance (cm), and rie and uh are the electron 
and hydrogen densities (cm"'') respectively, and these values are 
scaled for a circular area of sky of 20' radius (1257 arcmin ). The 
range for each norm found by looking at all 13 Suzaku background 
regions and using a column density range between 0.32-0.52 x 10^^ 
cm"^ was 0.01-0.02 cm"'^ for the O.lkeV norm and 0.0009-0.0017 
cm~^ for the 0.6 keV nonn. These ranges include the effect of ac- 
counting for stray light emission in the background regions which 
is modelled using the ray tracing simulator XISSIM and subtracted. 

When performing the initial spectral fits to the cluster emis- 
sion we use the best fitting galactic APEC norm values for the 
galactic background obtained using the average column density of 
0.42x lO^'^ cm~'^, and then later investigate the systematic error on 
the fits of varying all of the background components (including the 
CXB and the NXB) by their quantified uncertainties while letting 
the column density be a free parameter. This is described in detail 
in section [431 



3.5 Stray Light and PSF spreading 

Suzaku observations suffer from stray light effects due to X-rays 
passing through the telescopes along undesired optical paths, such 
as reflection purely from the secondary mirror and reflection from 
the back surface of the primary mirror dSerlemitsos et al.ll2007h . 
This results in X-rays being detected from bright sources outside 



heasarc.gsfc.nasa.gov/docs/rosat/ros_xselect_guide/xseIect_ftooIs.htmI 



of the field of view of the telescope. In addition there is also light 
spread from adjacent regions due to Suzaku's broad PSF. 

We estimate the relative contributions into each annulus due 
to the PSF and stray light as follows. Using the best fit Suzaku sur- 
face brightness profile found later in section |6]for each annulus as 
input to the ray tracing simulator XISSIM using the same number 
of photons (2x10''), we measure the number of photons scattered 
into the other annuli, and calculate the fractional contamination in 
each annulus from all of the other annuli. This is shown in table 
[3] An error in the calculation of the mixing fractions in G09 led 
to the contributions entering the outer annuli from the inner annuli 
being underestimated in table 2 of G09, (for instance the contribu- 
tion from the core into the outermost annulus in G09 is incorrectly 
estimated to be 0.1 percent). The fractional contributions of stray 
light entering the outer annuli we calculate ar e similar to levels 
that have been measured in other Suzaku studies I Bautz et alj2009l : 



iHoshino et all20ld : lAkamatsu et alj201 ll : IWalker et alj|2012l) . 

We see that the majority of the emission observed in each an- 
nulus originates from that annulus, and the majority of the contam- 
ination originates from the neighbouring annuli which will have 
similar temperatures and so the effect of this should be low. Con- 
tamination from the bright central 6' is non-negligible, and so when 
performing fits to the regions outside 6.0' we investigate the sys- 
tematic error of including in the background model a model of the 
emission reaching those regions from the central 6' (see section 
14.2b . This was achieved by using XISSIM to simulate the emission 
in all of the pointings using the spectrum, image and flux of the 
central 6' as input. 



4 SPECTRAL ANALYSIS 

To investigate the cluster emission we divide the pointings into the 
annuli shown in Fig.[T] with the annuli between 0.0'-2.5', 2.5'-6.0', 
6.0'-8.8', 8.8'-11.7', 11.7'-14.5', 14.5'-17.3', 17.3'-20.2', 20.2'- 
23', all of which are at least as wide as the Suzaku PSF (2.5') to 
ensure that the majority of the flux measured in each annulus orig- 
inated from that annulus. We have already concluded earlier that 
there is no statistically significant emission above the background 
level for radii greater than 23'. 



4.1 Azimuthally averaged deprojected fits 

We initially fit each direction in the cross separately but found no 
statistically significant azimuthal variation in the temperature pro- 
files, so we focus on azimuthally averaged profiles. We find the 
azimuthally averaged deprojected temperature and density profiles 
over the available pointings. We perform fits in the 0.7-7.0 keV 
band and allow the metallicity to be free for the central 3 annuli, 
but fix it to 0.3 Z0 for the outer 5 annuli which is commonly found 
and was used in G09 and Ella. For the central 3 annuli we find 
metallicities of 0.37l^ EJ^ Zq, 0.34l^-^g Zr and0.2+ ^ i Zq, con- 
sistent with the values found in lSnowden et al.l ^20081) using XMM- 
Newton. 

We allow the column density to be a free parameter for each 
region to accommodate for its spatial variation. The values for the 
column density found by letting this parameter be free are shown 
in Fig.|4] and are all within the range indicated by the LAB survey. 
The spectral fits are shown in Fig. lAll 

The deprojected temperature, density and entropy 
(S=kT/ne'''') profiles are shown in Fig. |5] We model each an- 
nulus as the sum of the absorbed APEC components describing 
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Table 3. Percentage contribution of flux in the rows' annulus from the columns' annulus due to PSF spreading and stray light. 



0.0'-2.5' 2.5'-6.0' 6.()'-8.8' 8.8'-11.2' 11.2'-14.5' 14.5'-17.3' 17.2'-20.1' 20.1'-23.0' 



0.0'-2.5' 


96.6 


3.4 


0.01 


0.002 


0.0003 


0.0002 


0.0002 


5x10^5 


2.5'-6.0' 


38.5 


60.5 


1.0 


0.03 


0.008 


0.003 


0.002 


0.001 


6.0'-8.8' 


16.8 


14.1 


64.6 


4.0 


0.43 


0.1 


0.04 


0.02 


8.8'-11.2' 


12.7 


7.0 


8.9 


64.9 


5.8 


0.6 


0.1 


0.03 


11.2'-14.5' 


5.5 


4.4 


1.8 


12.9 


64.7 


10.1 


0.5 


0.1 


14.5'-17.3' 


1.0 


1.4 


1.0 


1.7 


14.7 


72.1 


7.6 


0.4 


17.2'-20.I' 


2.9 


1.1 


0.6 


1.2 


2.1 


19.9 


67.0 


5.2 


20.1'-23.0' 


8.7 


1.9 


0.4 


0.7 


1.4 


2.5 


16.1 


68.4 



emission originating from the annulus itself and from the shells 
outside it, scaling each APEC component according to the volumes 
projected onto the annulus from the external shells assuming spher- 
ical symmetry. This emulates the XSPEC mixing model PROJCT, 
and allows data from different detectors to be fitted simultaneo usly 
(this is the same m ethod described in lHumDhrevetai]|201ll and 
IWalker et al.ll2012h . We rebin the temperature profile with more 
course binning owing to the limited spectral quality of the data, 
which is achieved by tying together the temperatures of the two 
annuli between 8.8' and 14.5', and the two annuli between 17.3' 
and 23.0'. We fit all 8 annuli simultaneously to propagate errors. 

4.2 Systematic errors 

We now turn to investigate systematic errors in the measurements. 
We need to quantify the systematic error on the temperature and 
density profiles resulting from our uncertainty in the background 
model used. As we found in section [3T1 the uncertainties in the 
variation of the column density increases the uncertainty in the 
level of the galactic emission, and we found that the normalisations 
of the 0. 1 keV and 0.6 keV APEC components could only be con- 
fined within the ranges of norm(0. lkeV)=0.010-0.020 cm~^ and 
norm(0.6keV)=0.0009-0.0017 cm"^ The CXB is also expected to 
vary spatially owing to variations in the number of unresolved point 
sources in each annulus, and the one sigma variations from the ex- 
pected 2-10 keV band flux of 1.87±0.13x 10"" erg cm"^ s"^ 
deg~^ for each annulus have already been calculated in table |2] In 
addition to this th e NXB level has an uncertainty of ±3 percent 
iTawa et alj j2008l) . 

Previous studies of cluster outskirts (for instance G09 and 
lAkamatsu et alj I2OI ih have quantified the effect of these back- 
ground uncertainties by varying each parameter by its one sigma 
error and measuring the effect on the temperature and densities. 
This approach is overly simplistic and underestimates the true sys- 
tematic eiTor caused by the uncertainty on the background model. A 
more realistic approach requires that all of the background compo- 
nents are varied simultaneously through their expected variances, 
and the effect on the temperature and density profiles measured. 

To achieve this we follow the method of Walk er et al.l ( [2012h 
and produce 10000 realisations of the background model (the 
galactic APEC norms, the CXB norms and the NXB level) dis- 
tributed by their expected variations which we have already found, 
and accounting for covariance between the background parameters. 
We then perform the deprojected fit for each realisation of the back- 
ground, thus folding our uncertainty in the background modelling 
through the deprojection to completely propagate the errors. From 
the resulting 10000 temperature and density profiles we obtain, the 
one sigma systematic errors are found for each annulus by finding 



the ranges which encompass 68 percent of the values around the 
peak value. These systematic errors are plotted as the solid black 
lines in Fig.|5] 

We then perform the deprojected fits including the stray light 
emission from the central 6' (the bright core) as a background com- 
ponent, and plot the effect of this on the temperature and density 
profiles as the green line in Fig.|5] We find that taking into account 
the stray light is a small effect which is well within the statistical 
errors of the fits. We vary the cont amination on the optical blocking 
filter (OBF) by ±10 percent (as in lAkamatsu et al.l201ll) . using the 
ftool XISCONTAMICALC to modify the ARFs, and find that this has 
a negligible effect which is much smaller than the statistical eiTors 
on the temperatures and densities. 

4.3 Discussion of temperature, density and entropy profiles 

As mentioned in Ella, the temperature profile found in G09 
is significantly lower than tha t obtained using X MM-Newton 
jSnowden etal]r2008) and Swift jMoretti et al1l201 ih outside 6 ar- 
cmins, leading to a lower mass and r2oo measurement than has been 
obtained with XMM-Newton and Swift. We find that our tempera- 
ture profile in the range 6.0'- 14.5' is significantly higher than that 
obtained in G09 and is consistent with the Swift and XMM-Newton 
results. This is shown in the top panel of Fig. |5] This is a result 
of a combination of more accurate background modelling than that 
used in G09, and also of letting the column density be a free pa- 
rameter. G09 fixed the column density to 0.37 x 10^^ cm~^, which 
was found to be the best fitting value when the column density was 
allowed to be free for the the central region ( within 2.5'). 

For comparison, the Swift results of iMoretti et al.l ( 1201 ih 
fix the column density to the mean LAB survey value of 
0. 418x10^^ cm~^ for t he entire cluster. The column density used 
in lSnowden et al. U2008I) for the XMM-Newton results was allowed 
to be a free parameter for each annulus studied. 

The deprojected density profile agrees well with that obtained 
with ROSAT in El la out to 17', and is consistent with the 90 per- 
cent upper limit beyond this when the systematic errors in the back- 
ground modelling are taken into account. As described in detail 
later in section |6] the ROSAT surface brightness (and therefore 
density) in Ella appears to have been underestimated in the out- 
skirts due to inconsistent point source removal which has caused 
the background to be overestimated in Ella. As no emission is 
found outside 23' the excess cannot be the result of projected emis- 
sion from outside the 20.l'-23' annulus. Including the emission 
spread into the annulus from neighbouring annuli due to the PSF 
has no statistically significant effect. 

2/3 

We find the entropy (S=kT/ne ) profile using the depro- 
jected temperatures and deprojected densities as shown in the bot- 
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Figure 5. Top: Deprojected, azimuthally averaged temperatur e profile from 
Suzak u shown with black points. Blue points ar e from Swift jMoretti et al] 
1201 ih . the red points are from XMM-Newton (Snowden et al. 2008) and the 
pink points from from BeppoSAX (De Grandi & Molendi 2002), showing 
good agreement with previous observations. The projected temperatures ob- 
tained in G09 in the 6.0'-13.5' region are shown as the cyan points, and are 
significantly lower as described in Ella. Middle: Deprojected density pro- 
file from Suzaku is in black. The deprojected ROSAT density profile from 
El la is in red, where the outer 2 density values are 90 percent upper limits. 
Bottom: Enfi'opy profile obtained using deprojected densities and depro- 
jected temperatures. The entropies corrected by the clumping factors found 
later in Fig. [Tjare shown as the blue points. The dashed black line shows 
the theoretical p rediction (S cx: r^'^) from the numerical simulations of 
IVoit et al"] j2005). One sigma systematic errors owing to background mod- 
elling are shown as the solid black lines for the temperature, density and 
entropy profiles. The solid green line show the best fits obtained when the 
stray light from the central 6' is subtracted from the cluster spectra. 



torn panel of Fig. [5] and co mpare it t o the t heoretical predic- 
tion of S oc r^^ obtained in IVoit et alj ( l2005h assuming purely 
gravitational hierarchical structure formation (and neglecting non- 
gravitational processes such as cooling and feedback). The higher 
temperature found in the 6.0'- 14.5' region compared to G09 causes 
the entropy profile to be consistent with the S oc r^'^ relation 
out to greater radius than in G09 (15' compared to ~10'). Out- 
side 15' the entropy profile flattens and is much lower than the 
S oc r^'^ relation, which is due to the low outer temperature (~ 
2keV). This entropy profile flattening has been observed in other 
cluster s observed to the virial radius with Su zaku, namely Abe ll 
1689 jKawaharada et al."2010^, Abell 1795 jSautz et alj |2009|). 



the Perseus clus ter (Simionesc u et al. (201 lall) and Abell 2142 
jAkamatsu et alj|201 1) ), and may indicate that the ICM is out of 
hydrostatic equilibrium due to the continual infall of matter onto 
the cluster (which has not been virialised yet) as it continues to 
form. In section |5] we find that r2oo (the boundary within which 
we expect the cluster to be virialised) lies between 16.7' and 18.4', 
which is consistent with the suggestion that the cold, low entropy 
gas we find outside 17' is out of hydrostatic equilibrium and has 
not yet been virialised as it accretes onto the cluster 

In Fig. [6] we compare the entropy profile (scaled by the 
average ICM temperature) with the entro py profiles found for 
other clusters investigated by Suzaku to r;>nn. lSato et al. Il l2012h has 
shown that when scaled in this way the entropy profiles are con- 
sistent with a low level of scatter, but that the PKS 0745-191 re- 
sults from G09 are lower and inconsistent with the trend. Our new 
entropy profile shows much stronger agreement with the general 
trend, resulting mostly from the higher temperatures we find be- 
tween 6.0'- 14.5'. 



5 MASS ANALYSIS 

We assume that the ICM is spheri cally symmetric and is in hydro- 
static equilibrium, so that it obeys Jvikhlinin et al.l2006l) . 

M« r) = -3.68 X 10"Mer(Or ( + (3) 

\ a m r a In r / 

Furthermore, we assume that the total mass density obeys a 
Navarro Frenk White (NFW) profile (Navarro et al. 1997) (this as- 
sumption is strongly motivated by numerical simulations of hierar- 
chical structure formation) 

Po 



p(r) 



po 



(4) 



r/rs{l + {r/r^)y 
200/9cC^oo/3(«n(l + caoo) - C2oo/(l + C200)) (5) 



where pc is the critical density of the universe, and 
C20o=r2oo/i's- An i n G09, we follow a mass analysis technique sim- 
ilar to that used in Schmidt & AUeril 1I2OO7V Using the deprojected 
density profile and the NFW mass profile we calculate the tem- 
peratures in each annulus, beginning with the outermost annulus, 
and compare them to the observed temperatures. The results are 
unchanged if we start with the innermost annulus. The best fit 
NFW mass profile corresponds to that which produces tempera- 
tures which match the observed temperatures most closely, so that 
it minimises the statistic defined as 

(^Tcalculated^i Tactual, 



(6) 



We generate 10000 realisations of the density profile dis- 
tributed by the combined statistical and systematic errors shown 
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Figure 6. Azimuthally averaged entropy profiles scaled by the average ICM temperature for clusters investigated by Suzaku to r200 . showing a low level of 
scatter The black crosses are our new results for PKS 0745-191. The rad ial error bars for the other clusters are omitted to improve figure clarity. Green is for 
Abell 2029 (Walker etal. 2ol l), blue is for Abell 2 142 jAkamatsu et ali 201 1). red isforAbell 1689 (Kawaharada et al. 2010), cyan is for Hydra A jSato et alj 
|20T^ . grey is for Abell 1413(Hoshino et al.l2O10h and pink is for Perseus (NW direction only from lSimionescu et al. (20lla1l) . The solid black line shows the 
S oc r^ i relation from iVoit et al..a005,) . 



in Fig. |5] and repeat the mass analysis for each profile, calculat- 
ing the cumulative gas mass fraction profile for each iteration. The 
resulting gas mass fraction profile is shown in Fig. |7] along with 
its 1 G errors. To resolve the central annulus, which is needed for 
constraining the NFW profile, we use the S wift temperatures an d 
densities for the regions within 2' obtained in lMoretti et al.l ( l201lh . 

We find that our higher Suzaku temperatures and more accu- 
rate eiTor measurements compared to G09 allow a much better fit 
to an NFW profile to be obtained, with a best fit x^/d.o.f of 9.0/7 
(compared to 39/7 obtained in G09) when we used the entire tem- 
perature profile except for the outermost temperature annulus. We 
do not include the outermost temperature annulus in the mass anal- 
ysis because when we perform a simple mass analysis using linear 
gradients between the temperature points in the hydrostatic equi- 
librium equation we find that the hydrostatic mass unphysically de- 
creases in the outermost annulus due to the low temperature, indi- 
cating that the ICM is out of hydrostatic equilibrium there. Includ- 
ing this low temperature outer annulus gives a much poorer fit to 
an NFW profile (x^/d.o.f = 14.0/8), which, together with the flat- 
tening of the entropy profile there, adds further evidence that the 
ICM is out of hydrostatic equilibrium there. This is similar to the 
case of Abell 1689 studied in Kawaharada et al. ( 2010), where the 
directions in which the entropy profile flattened were found to have 
mass profiles which unphysically decreased in the outskirts, indi- 
cating the breakdown of hydrostatic equilibrium. 

When the 14.5'-17.3' temperature annulus was removed from 
the mass analysis the NFW parameters were not significantly 
changed, resulting in the same 1200 and M200 values. This suggests 



that in the I4.5'-17.3' annulus the ICM is still in hydrostatic equi- 
librium, and that it is only outside ~17.3 arcmins that deviations 
from hydrostatic equilibrium or non- thermal pressure support be- 
come significant. 

We find M2oo=9.6ll:|^ xlO" Mq, C2oo=5.3lg:^, 
r,=370l5o° kpc, and r20o=2.0to\ Mpc (17.4tJ? arcmins), 
which are consistent with t he res ul ts obtained for PKS 0745- 
191 in Pointecoutea u et al.1 i2005h . IVoigt & Fabiaiil ('2006) and 
Schmidt & Allen (20071) . and with the scaling relations of 
Amaud et al. (j2003), all of which are tabulated in table |4] It must 
be stressed however that the Chandra results of I Voigt & FabiarJ 
(f2006) and I Schmidt & AUenl ( [2007 ) used temperature profiles 
extending out to o nly ^400 kpc (~3.5 arcm ins), while the XMM- 
Newton results of IPointecouteau et al.l ( 12005.) used a temperature 
profile which only extended to ~1000kpc (~9 arcmins), and that 
these were extrapolated outwards to find M200. 

We see therefore that the lower temperatures obtained in G09 
due to incorrect background modelling caused the mass determina- 
tion in G09 to be too low and the NFW concentration parameter to 
be unusually higher than previous measurements. 

We find that the gas mass fraction rises with increasing radius 
and rises ab ove the mean cosmic baryon fraction obtained using 
CMB data in lKomatsu et alj t20I ll ) from around r2oo onwards, but 
is lower than (but still consistent with) the cumulative gas mass 
fraction profile for the Perseus cluster in lSimionescu et al. (201 Ia)l 
(scaled to the r2oo value we find), as shown in Fig. |7] . If we as- 
sume that the cumulative gas mass fraction is overestimated due 
to gas clumping alone then the clumping factor (by which the gas 
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Figure 7. Top panel: Gas mass fraction profile; the horizontal dashed black line shows th e cosmic mean baryon fraction determined in lKomatsu et alj fcOllI) . 
The blue points are those obtained for the Perseus cluster in ISimionescu et al. (201 lail scaled to the r2oo value we find for PKS 0745-191. Middle panel: 
Clumping factor by which the gas density must be overesti mated if the true cumulativ e gas mass fraction does not increase above the mean cosmic baryon 
fraction. The blue fine is the clumping factor detennined in lSimionescu et al. (201 lal for the Perseus cluster. Bottom panel: The black line shows the factor 
by which the total mass needs to be underestimated for the cumulative gas mass fraction to remain below the mean cosmic baryon fraction, with the red lines 
sho wing the one sigm a error. The green and blue sold lines show the expected mass underestimates from the the numerical simulations of lLau et"al] j20O9l) 
and lBums et alj j2010l) respectivelv. while the dashed green and blue lines show the upper limit from these simulations. All panels: The best fit is in black and 
the red lines show one sigma eiTors. r2oo is shown by the vertical dashed black line and its one sigma errors by the vertical dashed red lines. 



density is overestimated) which is required for the cumulative gas 
mass fraction not to exceed the cosmic mean baryon fraction is 
s hown in the middle panel o f Fig. [7] and is lower than that found 
in lSimionescu et al. (201 lail for the Perseus cluster, but is still we ll 
within range indicated by the simulations o f lNagai& LaiilfeOllh . 

Assuming that the gas density has been overestimated by the 
clumping factors calculated in Fig. [7] we recalculate the entropy 
profile in the outskirts as the blue points in the bottom panel of Fig. 
[3] The clumpi ng corrected entr opy profile agrees with the power- 
law relati on of IVoit et al.l ( 2005h out to our estimate of r2oo (as was 
found in ISimionescu et al. (201 la"ih . but outside r2oo the entropy 
profile still flattens and lies below the powerlaw prediction, sug- 
gesting that outside the virial radius the entropy profile flattening 
must be caused by something in addition to clumping. The most 
likely additional cause for the entropy profile flattening outside r2oo 
is that the ICM is out of hydrostatic equilibrium, and that we are 



seeing cold gas which is still accreting onto the cluster which has 
yet to be virialised. 

Following the analysis of ISimionescu et al. (201 lal we then 
compare the pressure profile with the universal pressure profile 
found in lAmaud et all ( I2OI0I) by fitting to clusters explored with 
XMM-Newton out to ~ rsoo and using simulations beyond this, and 
see the effect that correcting for c lumping has (Fig. |8). As was 
found in ^Simionescu et al. (201 lal the pressures in the outskirts 
are systematically higher than the extrapolated universal pressure 
profile from around rgoo outwards, but they are still consistent with 
the universal pressure p rofile within the range obtained from hy- 
drodynamic simulations jLapi et al.l20T2h . shown by the cyan lines 
in Fig. [8] Taking into account the clumping factor reduces the pres- 
sures within r2oo so that they become m ore consistent with the uni- 
versal pressure profile, as was found in lSimionescu et al. (201 lal . 
but outside r2oo correcting for clumping causes the pressure to lie 
below the universal pressure profile. This adds further evidence to 
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the hypothesis that outside 1200 the gas is out of hydrostatic equi- 
librium, and that we are seeing cold gas which is accreting onto the 
cluster. 

It is also possible that the gas mass fraction is being over- 
estimated because we are underestimating the actual mass by 
using only the thermal pressure in the equatio n of hydrostatic 
equilibrium. Num erical simulations (for example iLau et all 1 20091 
iBums et^boiol) indicate that the pressure contribution from ran- 
dom gas motions increases with radius, leading to an underestimate 
of the true cluster mass if only the thermal pressure is used in the 
equation of hydrostatic equilibrium. In the bottom panel of Fig. [7] 
we show the radial profile of the factor by which the mass would 
need to be underestimated for the cumulative gas fraction not to ex- 
ceed the mean cosmic baryon f raction, and comp are t his with the 
simulated underestimates from iLau et aL I ( l2009t) and iBums et^ 
( l2010h . At r2oo the underestimate needed is ~15 percent, which 
is highly consistent with the simulated underestimates. Outside the 
r2oo the underestimate rises but remains consistent with the simu- 
lations of iLau et al.l h00% . 

We therefore conclude that the high gas mass fraction can be 
explained either by gas clumping or by an understimate of the to- 
tal mass which is within the expected range from numerical sim- 
ulations. In reality it is likely that both factors are contributing to- 
gether, but it is not possible to determine the relative contributions 
of each individually. 



6 COMPARING ROSAT AND SUZAKU SURFACE 
BRIGHTNESS PROFILES. 

EI la took the projected densities obtained in G09 and used these to 
calculate a MEKAL norm which they then folded through a ROSAT 
PSPC response to obtain the projected surface brightness profile 
that the Siizaku results of G09 predict would be obtained using 
ROSAT PSPC. When the surface brightnesses were compared in 
EI la a significant discrepancy was found outside 14 arcmins, with 
the Suzaku surface brightness being higher than the ROSAT at a 
level of 7.7a. 

We repeat the ROSAT imaging analysis of Ella (us ing the 
Extended Source Analysis Software o f ISnowden et al.ll994 ') but in- 
stead use a constant threshold count rate for resolving point sources 
of 0.003 counts per second in the R47 ba nd when using the ftool 
DETECT (as used in lEckert et al. (201 Ibih . meaning that the CXB 
is resolved to the same extent across the entire detector The PSF of 
the ROSAT PSPC increases strongly away from the on axis posi- 
tion, meaning that the point source resolving ability is not uniform 
across the detector, and that it is able to detect fainter point sources 
nearer to the on-axis position. In order for the CXB to be resolved 
to the same level uniformly across the entire detector it is necessary 
to resolve the point sources down to the same threshold flux across 
the entire detector. 

From Fig. 1. in El la, which shows the point sources removed 
in their analysis, it seems from the greater number of resolved point 
sources closer to the on axis position that a uniform threshold point 
source flux has not been used. The exact threshold used in remov- 
ing the point sources is not mentioned in EI la, but we find by re- 
peating their analysis that the mask used can be produced using 
the ftool DETECT with a detection threshold of 3 a, and that us- 
ing this reproduced the findings of El la. The masked PSPC image 
is shown in Fig. IA3I along with the circle at the r2oo value deter- 
mined in G09, which shows the same point source masking shown 
in Fig. I of Ella. The use of a constant significance means that 



more point sources have been removed in the cluster emission re- 
gion (nearer the on-axis position) than in the background region 
used in El la, meaning that the CXB has been resolved more in the 
cluster emission regions. To make this clear, the point sources in 
Fig. IA3I are colour coded according to their count rates, showing 
that the point sources resolved in the cluster region are fainter than 
can be resolved in the background region. The background level 
used in Ella is therefore an overestimate. 

When we use a constant threshold count rate of 0.003 counts 
per second in the R47 band to resolve point sources, the resulting 
surface brightness profile in the outskirts, extracted in the same re- 
gions as our Suzaku observations, is shown as the blue points in 
Fig.|9] We then use our best fit projected Suzaku results to produce 
the projected surface brightness profile we would expect to see with 
ROSAT in the R47 band. This is shown as the black points in Fig.|9] 
where the statistical error on the fits and the systematic error in the 
background modelling have been added to produce the error bars 
shown. We find complete agreement between the ROSAT results 
and the Suzaku results. 

The dominant cause of the discrepancy between the Suzaku 
results of G09 and the ROSAT results of EI la is the underestimate 
of the background level in G09, as when a constant flux threshold is 
used for the ROSAT analysis the predicted surface brightness from 
the G09 results is still higher than the ROSAT measurement at the 
5a level in the outskirts. 

The results in Ella are shown as the red points (which used 
the same binning as G09) in Fig.|9] and we see that the use of a con- 
stant significance in removing point sources (rather than a constant 
flux) in El la has caused the CXB to be resolved more near the cen- 
tre of the pointing than in the background regions, thus causing the 
background to be overestimated and underestimating the surface 
brightness in the outskirts of the cluster. This subsequently means 
that the densities found in the cluster outskirts in Ella (which are 
shown in Fig.|5]( are also likely to be underestimated. 

Ella used the density profile obtained with ROSAT and the 
temperature profiles obtained with XMM-Newton, BeppoSAX and 
Swift to perform a similar mass analysis to the one we present and 
their results are tabulated for comparison in table |4] Whereas the 
masses using the BeppoSAX temperatures and the ROSAT densi- 
ties are consistent with our results, those obtained using the XMM- 
Newton temperatures and Swift temperatures are higher than our 
results. The ROSAT mass analysis of El la used the density profile 
out to 17 arcmins, which from Fig.|5]is completely consistent with 
our Suzaku density profile, so the differences in the mass determi- 
nations are due to the differences in the temperature profiles, as is 
also described in Appendix A of El la. 



7 SUMMARY 

We have used new Suzaku observations of the background regions 
surrounding PKS 0745-191 to perform a thorough analysis of the 
error in the background modelling performed, allowing a more ac- 
curate measurement of the thermodynamic properties of the ICM 
out to 23', beyond which we do not detect any statistically signif- 
icant emission. We have highlighted and corrected the problems 
with the background modelling used in G09. 

We find a temperature profile which is consistent with those 
already obtained using XMM-Newton and Swift, and conclude that 
the low temperatures found in G09 in the 6.0'- 14.5' region were the 
result of incorrect background modelling. Our higher temperatures 
cause the entropy profile to obey the powerlaw relation, S oc r^ '^. 
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Figure 8. Solid line shows the universal pressure profile of lAmaud et aP told) out to rsoo, which is extrapol ated outwards as t he dashed hne. Cyan hnes shows 
the range in the universal pressure profile from hydrodynamical simulations of relaxed clusters presented in iLapietalJiioTi) . Black point show the pressures 
(kUeT) using our deprojected densities and deprojected temperatures, which agree with the universal pressure profile within the range of the simulations. Blue 
points show the pressures corrected by the gas clumping factors found in Fig. [T] which brings the pressures into better agreement with the universal pressure 
profile inside r2oo > but which causes them to lie below the universal pressure profile outside r2oo ■ 



Table 4. Comparing our mass determination with previous estimates. The resu lts from lPointecouteau et al.lj2005l)JVoigt & FabianH200^ and lSchmidt & Allenl 
were taken from the table compiled in IComerford & NataraiarJ bOOVh . 



C200 


M200 /IQi* Mq 


Reference 


Comments 


c 1+0.6 
-'■■^-0.9 


n + 1.9 
^■°~1.0 


This work 


Neglects outermost temperature measurement (x'^/d.o.f = 9.0/7) 


, ,7+0.40 
J-i^_0.40 
5 46+3-22 


lo.oiii 

^ y-|-52.2 

lO.O+i'j 

<: 4+0.6 


Pointecouteau et al. (2005) 
Voiat & Fabian (2006) 
Schmidt & Allen a007) 
Amaud et al. (2005) 
Georae et al. (2009) 


Uses XMM-Newton temperature profiles extending out to ~1000 kpc (~9 arcmins) 
Uses Chandra temperature profiles extending out to ~400 kpc (~3.5 arcmins) 
Uses Chandra temperature profiles extending out to ~400 kpc (~3.5 arcmins) 

-1-0 1 

M-T scaling relations (r200=1999_77kpc=17.6_Q Yarcmins) 
Underestimated background modelling (x^/d.o.f = 39/7) 


14^+0.25 

-"■^^-0.25 
4 BQ+0.89 
^■°^-0.89 
c-,+0.26 
-'-''^-0.26 


17.48ti;i 

12 03+30'' 
1Z.UJ_3 04 

18.0012-05 


Ella 
Ella 
Ella 


Uses ROSAT density profile and XMM temperature profile, (x^/d.o.f = 57/23) 
Uses ROSAT density profile and BeppoSAX temperature profile, {^ld.o.f = 27/18) 
Uses ROSAT density profile and Swift temperature profile, (x^/d.o.f = 61/19) 



to higher radius than in G09 before flattening and decreasing above 
15' (~ 1.7 Mpc). This brings the average temperature scaled en- 
tropy profile into strong agreement with those from other clusters 
investigated with Suzaku out to r2oo • 

The higher temperatures found in the 6.0'- 14.5' region im- 
prove the fitting statistics for modelling the ICM as being in hy- 
drostatic equilibrium in an NFW potential from x^/d.o.f =39/7 ob- 
tained in G09 to x^/d.o.f =9.0/7 using all of the temperature profile 
except the outer aimulus. The fit to an NFW profile is much poorer 
if the outermost temperature annulus is included, further indicating 
that beyond ~ 17' the ICM is out of hydrostatic equilibrium. The es- 
timates of M200, C200 and r2oo are consistent wi t h prev ious results 
and with the scaling relations of lAmaud et al.l (l2005h . as shown 
in table |4l thus resolving the problem with the G09 measurements 



which found unusually small M200 and r2oo values. We find that 
r2oo lies in the range 16.7-18.4 arcmins (1.9-2.1 Mpc). This means 
that we have probed the ICM of PKS 0745-191 to at least 1.25r2oo 
= 23', beyond which we find no statistically significant emission. 

We find that the gas mass fraction increases significantly 
above the mean cosmic baryon fraction from around r2oo onwards, 
and is consistent with the profile found in lSimionescu et al. (2011 ail 
for the Perseus cluster. This suggests that the ICM may be clumped 
in the outskirts, causing the gas density to be overestimated. Cor- 
recting for the clumping factor (if we assume that the true gas mass 
fraction does not increase above the mean cosmic baryon fraction 
purely because of gas clumping) causes the entropy profile to obey 
the powerla w relation S oc r^'^ out to our estimate of r2oo (as was 
observed in Isimionescu et al. (2011a)l for the Perseus cluster) but 
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Figure 9. Comparing the projected ROSAT PSPC surface brightness pro- 
file in the outskirts found using a constant point source thi'eshold count 
rate of 0.003 counts per second (blue points) with our projected Suzaku re- 
sults folded through the ROSAT PSPC response in the R47 band (shown in 
black). A small horizontal offset has been included between the points to 
increase the clarity. The eiTors on the Suzaku results add together the one 
sigma statistical errors of the fits and the one sigma systematic errors of the 
background modelling. The results from El la are shown in red, where the 
surface brightness in their outer bin is a 90 percent upper limit. 



above r2oo the entropy profile still flattens and lies below the pow- 
erlaw prediction, suggesting that something other than clumping is 
responsible for the entropy profile flattening there. The most likely 
cause is that outside r2oo the ICM is out of hydrostatic equilibrium 
and that we as seeing cold gas which has not yet been virialised 
as it accretes onto the cluster. This hypothesis is further supported 
by comparing the clumpi ng corrected pressure profile with the uni- 
versal pressure profile of lAmaud et al.l dlOlOl) . Between 1500 and 
r2oo the clumping correction reduces the pressure and improves 
the agreement with the un iversal pressure profile, as was found in 
ISimionescu et al. (201 lal However outside r2oo the clumping cor- 
rection causes the gas pressure to lie below the extrapolated univer- 
sal pressure profile, suggesting that outside r2oo we are seeing cold 
gas which is accreting onto the cluster. 

It is also possible that the observed gas mass fraction exceeds 
the cosmic mean baryon fraction in the outskirts because the total 
mass is being underestimated by using only the thermal pressure 
in the equation of hydrostatic equilibrium. The total mass underes- 
timate needed for the cumulative gas mass fraction to stay below 
the cosmic mean baryon fraction agrees with the expec ted under- 
estim ate of the total mass from numerical simulations jLau et all 
I2OO9I) . which is caused by the increasing contribution of non- 
thermal pressure in the outskirts which cannot be directly measured 
at present. In reality it is likely that a combination of gas clumping 
and non-thermal pressure support is responsible for the high gas 
mass fraction and the flattening of the entropy profile in the out- 
skirts. 

Our results indicate that a constant threshold flux was not 
used for resolving point sources in ROSAT PSPC analysis of El la. 



which has caused greater point source extraction nearer to the on- 
axis position due to the increase in the PSF (and the decrease in 
the point source detection sensitivity) away from the centre of the 
PSPC detector. This means that the CXB was not resolved uni- 
formly across the PSPC pointing, with the threshold flux for the 
point sources removed in the background regions being higher 
than it was for the cluster emission regions. This has caused the 
background to be overestimated, thus underestimating the surface 
brightness and density in the cluster outskirts. 

We predict from our Suzaku spectral fits a projected ROSAT 
surface brightness profile which is consistent with that obtained 
with ROSAT when the errors in the background modelling are taken 
into account and when a constant point source threshold is used for 
removing point sources in the ROSAT PSPC image. 

Using our background pointings for PKS 0745-191 (which 
lies 3 degrees above the galactic midplane), we find no sta- 
tistically significant evidence for the presence of the absorbed 
bremsstrahlung background component (due to unresolved dM 
stars) found in Masui et al.l ( 12009 ) in their study of a background 
pointing in the galactic midplane. This su pports the hypothes is that 
this background component observed in iMasui et all ( |2009|) orig- 
inates from unresolved dM stars whose number density decreases 
rapidly away from the galactic midplane (as shown in Fig. 4 of 



rapidly away trom 
iMasui et al.i2009l) . 
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Figure Al. Spectral fitting for the cluster emission. In reading order we show the azimuthally averaged spectra used to generate the results in Fig. \5\ the 
corresponding annuH are 0.0'-2.5', 2.5'-6.0', 6.0'-8.8', 8.8'-11.7', 11.7'-14.5', 14.5'-17.3', 17.3'-20.2', 20.2'-23'. The data from each detector and each 
editing mode were fitted simultaneously and are added here for display purposes only. The black lines through the points represent the best projected fits 
(background plus cluster emission) which are shown for simplicity, while the lower black line shows the cluster emission. The green line shows the X-ray 
background level, and the red points show the NXB level which is subtracted from the X-ray spectra. The fit quality is the same for the deprojected fits but we 
show the projected fits for simplicity. 
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Energy (keV) 

Figure A2. Spectral fitting to all of the Suzaku background regions between 23'-41' simultaneously. The spectrum is well fit using two APEC components of 
solar metallicity and zero redshift to model the galactic emission; an unabsorbed component at 0.1 keV modelUng the LHB and an absorbed component at 0.6 
keV. An absorbed powerlaw of index 1 .4 is used to model the CXB . The NXB level which is subtracted from the X-ray spectra is shown overplotted as the red 
points. 
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Figure A3. Background subtracted, exposure corrected and smoothed ROSAT PSPC image in tlie R47 band following the procedure of Ella. The point source 
removal is exactly the same as that used in El la (compare the inner part with Fig. 1 in El la), and we find that this is achieved using a constant significance 
of 3 cr in the ftool DETECT when removing point sources. Due to the increase in the PSF with off axis angle, this means that the CXB has been resolved more 
deeply closer to the on-axis position, leading to an overestimate of the background in El la. The colours of the point sources correspond to their count rates; 
red point sources have count rates less than 0.002 cts/s, the greens ones are between 0.002-0.003 cts/s, and the blue ones are >0.003 cts/s. The white circle at 
the r2oo value determined in G09 is the same as that shown in El la in their Fig. 1 and allows a direct comparison. The background region used in El la is the 
pink annulus between 25' and 50'. 
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